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1 Introduction 



One of the main goals of this article is to present a somewhat general framework 
for the quantization of classical observables on a cotangent bundle which are poly- 
nomials at most cubic in momenta. This approach enables us to investigate the 
quantization of classically Poisson-commuting observables, and hence to tackle the 
problem of quantum integrability for a reasonably large class of dynamical systems. 

What should actually be the definition of quantum integrability is a long stand- 
ing issue, see, e.g., [37j. The point of view espoused in this paper is the following. 
Start with a complete set of independent Poisson-commuting classical observables, 
and use some quantization rule to get a corresponding set of quantum observables; 
if these operators appear to be still in involution with respect to the commutator, 
the system will be called integrable at the quantum level. 

Our work can be considered as a sequel to earlier and pioneering contributions 
[m m 13 1221 Ell that provide worked examples of persistence of integrability from the 
classical to the quantum regime. The general approach we deal with in this paper 
helps us to highlight the general structure of quantum corrections and to show that 
the latter actually vanish in most, yet not all, interesting examples. 

Returning to the general issue of quantization, let us mention that our choice 
of quantization procedure, which we might call "minimal" , doesn't stem from first 
principles, e.g., from invariance or equivariance requirements involving some specific 
symmetry. Although this "minimal" quantization only applies to low degree polyno- 
mials on cotangent bundles, it has the virtue of leading automatically to the simplest 
symmetric operators that guarantee quantum integrability in many cases. In order 
to provide the explicit form of the quantization scheme, hence of the quantum correc- 
tions, we need a symmetric linear connection be given on the base of our cotangent 
bundle. In most examples where a (pseudo-)Riemannian metric is considered from 
the outset, this connection will be chosen as the Levi-Civita connection. 

To exemplify our construction, we consider a number of examples of classical 
integrable systems together with their quantization. For instance, our approach 
for dealing with quantum integrability in somewhat general terms allowed us to 
deduce the quantum integrability of the Hamiltonian flow for the generalized Kerr- 
Newman solution of the Einstein-Maxwell equations with a cosmological constant 
first discovered by Carter 01101 El- Also does our quantization scheme leads us 
to an independent proof of the quantum integrability for Stackel systems originally 
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due to Benenti, Chanu and Rastelli [HIS]- 

The paper is organised as follows. In Section |21 we gather the definitions of the 
Schouten bracket of symmetric contravariant tensor fields on configuration space, 
M. We make use of Souriau's procedure to present, in a manifestly gauge invariant 
fashion, the minimal coupling to an external electromagnetic field; this enables us 
to provide a geometric definition of the so called Schouten-Maxwell bracket. The 
related definitions of Killing and Killing-Maxwell tensors follow naturally and will 
be used throughout the rest of the paper. We recall the basics of classical integrable 
systems, with emphasis on the Stackel class. The main objective of the present 
Section is then to revisit some classic examples of integrable systems involving Killing 
tensors. Naturally starting with the Jacobi system on the ellipsoid, we prove, en 
passant, that it is locally of the Stackel type, even allowing for an extra harmonic 
potential. This extends previous work of Benenti |[3| related to the geodesic fiow of 
the ellipsoid. Similarly, we show that the Neumann system is also locally Stackel. 
A number of additional examples, not of Stackel type, e.g., the Di Pirro system, 
and the geodesic flow on various (pseudo-)Riemannian manifolds such as the Kerr- 
Newman-de Sitter solution and the Multi-Centre solution are also considered. 

We introduce, in Section El a specific "minimal" quantization scheme for ob- 
servables at most cubic in momenta on the cotangent bundle T*M of a smooth 
manifold M endowed with a symmetric connection V, extending a previous pro- 
posal pi]. This quantization mapping is shown to be equivariant with respect to 
the affine group of (M, V). The computation of the commutators of quantum ob- 
servables is then carried out and yields explicit expressions for quantum corrections. 
We also provide the detailed analysis of quantum integrability for a wide class of 
examples within the above list. 

The concluding section includes a discussion and brings together several remarks 
about the status of the "minimal" quantization that has been abstracted from the 
various examples dealt with in this paper. It also opens some prospects for future 
investigations related to quantum integrability in the spirit of this work. 

Acknowledgements: We are indebted to Daniel Bennequin for several very 
interesting remarks, and to Brandon Carter for fruitful correspondence. Special 
thanks are due to Valentin Ovsienko for a careful reading of the manuscript together 
with a number of most enlightening suggestions. 
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2 Classical integrable systems 



2.1 Killing tensors 

Let us start with the definition of the Schouten bracket of two polynomial functions 
on the cotangent bundle {T*M, uj = d^i/\dx^) of a smooth manifold M. Consider two 
such homogeneous polynomials P = P*^ "*'=(a;)^ji ■ ■ - ^ik Q — Q^^"'^^{^)^ii ■ ■ - iii 
of degree k and ^ respectively; we will identify these polynomials with the corre- 
sponding smooth symmetric contravariant tensor fields P" = P*i'-*''(x)(9ij ®---®di^ 
and Q« = Q'^-'^{x)di^ ® ■ ■ ■ ® di^. 

The Schouten bracket [P'*,Q'*]5 of the two contravariant symmetric tensors P" 
and (of degree k and d respectively) is the symmetric contravariant (fc + £ — 1)- 
tensor corresponding to the Poisson bracket of P and Q, namely 

[P\Q% = {P,Q}K (2.1) 

Using the the Poisson bracket {P, Q} = d^.PdiQ—d^-QdiP, and (|2.1|) . we readily 
get the local expression of the Schouten bracket of P" and QK If the manifold M is 
endowed with a symmetric connection V, the latter can be written as^ 

If M is, in addition, equipped with a (pseudo-)Riemannian metric, g, we denote 

by 

H = hg''i.i, (2.3) 

the Hamiltonian function associated with this structure. The Hamiltonian flow 
associated with H is nothing but the geodesic flow on T*M. 

A symmetric contravariant tensor fleld P" of degree k satisfying {if, P} = is 
called a Killing (or Killing-Stackel) tensor; using now the Levi-Civita connection V 
in ()2.2p . this condition reads 

y{ipii-iu) = 0. (2.4) 

^In this article the round (resp. square) brackets will denote symmetrization (resp. skew- 
symmetrization) with the appropriate combinatorial factor. 
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2.2 Killing-Maxwell tensors 
2.2.1 Souriau's coupling 

In the presence of an electromagnetic field, F, Souriau has proposed to replace 
the canonical symplectic structure, u, of T*M by the twisted symplectic structure 
Up = d^i A dx^ + ^Fijdx^ A dx^ . The (gauge-invariant) Poisson bivector now reads 

ttf = % Adi- \Fij % A d^^ . 

The Poisson bracket of two observables P, Q oiT*M is now 

{P, Q]f = npidP, dQ) = d^^Pd.Q - d^^Qd.P - P,, %P A d^^Q, (2.5) 

and the Schouten-Maxwell bracket of two polynomials P and Q is then defined by 

[P\Q%,F = {P,Q}1 

If the manifold M is endowed with a symmetric connection V, the Schouten- 
Maxwell bracket takes on the following form 

[P», Q^s,F = [P», Q^rr'^'-'d.. ® • ■ ■ ® d...... (2 6) 

with the expression ()2.2|) of the Schouten bracket [■ , ■]s. 

Suppose now that the manifold M is endowed with a metric g; the Hamiltonian 
vector field on {T*M,ijJf) for the Hamiltonian H given by ()2.3|) yields the the Lorentz 
equations of motions for a charged test particle moving on (M, g) under the influence 
of an external electromagnetic field P. 

A symmetric contravariant tensor field P" of degree k on (M, g) is now called 
a Killing-Maxwell tensor if {if, P}f = 0. The Killing- Maxwell equations then read, 
using ()2.6p . 

y{jpil...jfc) _ g ^ pi{il...ik-l p^k) _ g (2.7) 

where P/ = g^'^Fmi, in accordance with previous results [llj obtained with a slightly 
different standpoint. 

The conditions ()2.7p are of special importance for proving the classical and 
quantum integrability of the equations of motion of a charged test particle in the 
generalized Kerr-Newman background. 
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2.2.2 Standard electromagnetic coupling 

A more traditional, tliough equivalent, means to deal with the coupling to an elec- 
tromagnetic field, F = dA (locally), is to keep the canonical 1-form, a = C,idx^, on 
T*M unchanged, and hence to work with the original Poisson bracket {■,■}, but to 
replace the Hamiltonian ()2.3|) by 

#=|g^^(e.-A,)(e,-A,) (2.8) 

where the tilde makes it clear that the expressions to consider are actually polyno- 
mials in the variables C,i — Ai, for i = 1, . . . ,n; for example, if P = P^^"'^''C,ii • • • ^1^., 
then 

P = P^^-^^{^,^-AJ...{^,,-AJ. (2.9) 

The equations of motion given by the Hamiltonian vector field for the Hamil- 
tonian (j2.8j) on (T*M, da) are, again, the Lorentz equations of motion. 

The Schouten-Maxwell brackets and Schouten brackets for the electromagnetic 
coupling are related as follows via the corresponding Poisson brackets, viz 

{P,Q}f = {P,Q}. 

In this framework, a Killing-Maxwell tensor, P^, of degree k on (M, g) is defined 
by the equation {H,P} = 0. The resulting constraints are, again, given by ()2.7|) . 

From now on, and in order to simplify the notation, we will omit the jj-superscript 
and use the same symbol for symmetric contravariant tensors and the corresponding 
polynomial functions on T*M. 

2.3 General definition of classical integrability 

Let us recall that a dynamical system {/iA,uj,H) is (Liouville) integrable if there 
exist n = ^ dim independent Poisson-commuting functions Pi, . . . , P„ G C°°{M.) 
— that is dPi A ■ ■ ■ A dPn 7^ and {P^, Pf} = for all k,i = 1, . . . ,n — such that 
P, = H. 

We will, in the sequel, confine considerations to the case of cotangent bundles, 
{Ai = T*M, uo = dO) where 9 is the canonical 1-form, and of polynomial functions. 
Pi, ... , Pn, on T*M, that is to the case of n Schouten-commuting Killing tensors. 
Moreover, all examples that we will consider will be given by polynomials of degree 
two or three. 
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2.4 The Stackel systems 

These systems on {T*M,uj = dC,i A dx^) are governed by the Hamiltonians 

n 

H = Y,ci\x)^iC + M^')) (2.10) 

i=l 

where the i-th function /j depends on the coordinate only, and the functions a* 
are defined as follows. Let B denote a GL{n, R)-valued function defined on M and 
such that 

B{x) = (5i(x^)52(x2)...5„(x")) 

where the i-th column Bi{x^) depends on only {i = 1, . . . ,n); such a matrix will 
be called a Stackel matrix. Then take 

/a\x)\ 
a{x) = : 

\a"(x)/ 

to be the first column Ai{x) of the matrix A{x) = B{x)~^. 

The integrability of such a system follows from the existence of n quadratic 
polynomials 

n 

I, = J2M^){hC + M^l), i=l,...,n, H = h. (2.11) 

i=l 

We call Stackel potential every function of the form 

n 

U,{x) = J2 A\{^)M^'), ^ = 1, . . . , n; (2.12) 

i=l 

the potential appearing in the Hamiltonian is just Ui. 

One can check (see, e.g., pSJ, p. 101) that the n independent quantities Ii are 
such that 

n 

{/,, /^} = J2 i^e dsAl - d,A\)S,. {\et +ft), iy^m. 

s,t=l 

The relation A = B~^, gives the useful identity^ 

d,A] = -ClA^, Cl = j2^s^^ (2-13) 

s=l 

^The Einstein summation convention is not used. 
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which imphes 

A^dsAl^- A'^dsAl = 0, i ^ m, s,t = l,...,n (2.14) 
and therefore the so defined Stackel systems are classically integrable. 

Remark 2.1. Let us mention an interesting result due to Pars (see |2HI; P- 102): 
for a system whose Hamiltonian is of the form ()2.1()|1 . the Hamilton- Jacobi equation 
is separable if and only if this system is Stackel. 

Although these systems constitute quite a large class of integrable systems, they 
do not exhaust the full class. A simple example of a non-Stackel integrable system 
was produced by Di Pirro (see Section IT^ . 

2.5 The Jacobi integrable system on the elhpsoid 

Let £ C R""*"^ be the n-dimensional ellipsoid defined by the equation Qo{y,y) = 1 
where we define, for y, z & R""*'"'^, 

n 

Qx{y,z) = J2-^, (2.15) 

with < Co < ai < . . . < a„; the equations Q\{y, y) = 1 define a family of confocal 
quadrics. 

It has been proved by Jacobi (in the case n = 2) that the differential equations 
governing the geodetic motions on the ellipsoid, S, form an integrable system. The 
same remains true if a quadratic potential is admitted (see |2ZI)- The Hamiltonian 
of the system, prior to reduction, reads 

-t n n 

H{p,y) = ^Y.pl + ^Y.y" (2-16) 

where p,y & R""^^ and a is some real parameter. 

Moser has shown |2^ that the following polynomial functions 

Fa{p,y)=pl + ayl + J2 ^P-^'~^'^''^' with a = 0,l,...,n, (2.17) 

are in involution on (T*R""'"^, ^"=q dpa A dya). Those will generate the commuting 
first integrals of the Jacobi dynamical system on the cotangent bundle T*S of the 
ellipsoid. 
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Our goal is to deduce from the knowledge of ()2.17|1 the independent quantities 
in involution /i, ...,/„ on {T*£, d^i A dx') from the symplectic embedding 

given by Zi{p, y) = Qo{y,y) - 1 = and Z2{p, y) = Qo{p, y) = 0. 



Proposition 2.2. The restrictions F^lj^^^ = F^o l of the functions \2.11 ] Poisson- 
commute on T*S. 

Proof. We get, using Dirac brackets, 



{-Fa|T*£-5 -F/3lx'*£-} — {Fa, Fp}\j,, 



S 



1 (2.18) 

Fa}{Z2, F/s} — {Zi, Fp]{Z2, Fa]] 



for second-class constraints. Now, the denominator {^1,^2} = — 2 ^^^q (?/Q,/aQ,)^ 
doesn't vanish while {Zi,Fa} = A^paUal o.a)Zi — 4:{y'^/aa)Z2 is zero on T*S, for all 
a = 0, n. The fact that {Fa, F^} = completes the proof. □ 

The reduced Hamiltonian for the Jacobi system on the ellipsoid S is plainly 



0=0 a=0 



In order to provide explicit expressions for the function in involution Ji, . . . , J„, 
we resort to Jacobi ellipsoidal coordinates x^, . . . ,x"^ on S. Those are defined by 

Qx{y,y) = i-^y^ (2.20) 

where 

n n 

U,{X) = l[{X-x') and V {X) = 1[{X - aa) (2.21) 

1=1 a=0 

and are such that ao < a;^ < ai < < . . . < a;" < a„. The induced metric, 
g = Ylij^igij{x)dx^dx^, of the ellipsoid £ is given by 



5«J 



1 " 



X 
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and retains the form [22] 

n 



X 



i\2 



where gi(a;) 



1=1 



'T V{x^) 



(2.22) 



which is actually Riemannian because of the previous inequalities. We put for 
convenience g*(x) = l/gj(x). 

Using (j2.2Up and (j2.2H) . we find the local expressions ?/q,(x) via the formula 



n 



Qj/y iC 



2 i=l 



and then obtain the constrained coordinate functions 



2i/a / ^ 



given by the induced canonical 1-form Yl^=i = iYla=oPci^ya. 
The Hamiltonian ()2.19|) on {T*S, d^i A dx^) is then found to be 



1 " 

i=l 



.a=0 i=l 



X 



Note that the potential term is obtained from the large A behaviour 

Qx{y. y) ^ T^yl + yiY^ ^^yl + ■ ■ ■ 



a=0 



a=0 



which can be computed using relation ()2.2()|1 . One gets 



Q\{y,y) ~ J 



X 



a=0 



i=l 



+ 



(2.23) 



(2.24) 



(2.25) 



One relates the conserved quantities ()2.17p to their reduced expressions on T*S 
by computing, using (j2.24j) and (|2.23p . the expression of Fa\rp,^. One gets the 
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Proposition 2.3. The Moser conserved quantities {Fa\T*£)a=o n '^ctam the form 



F I 



where 

n n 

Gxii, x) = J2 9\^) n - ^')^' + a n (A - xO- (2.26) 

i=l jj^i i=l 

It is useful to introduce the notation crl{x) for the symmetric functions of order 
/c = 0,l,...,n — lof the variables (x^, . . . , x"), with the exclusion of index i, namely 

n 

HiX - = J2 (-l)'-^A"-V^i(x). (2.27) 

j^i k=l 

We note that, from the above definition, (Jq{x) = 1. 

It is also worthwhile to introduce other symmetric functions, ak{x), via 

n n 

n(A - xn = (-1)'A"-V,(x). (2.28) 

j=l k=0 

We thus have 

n 

Gx{^, x) = J2 (-1)^"'A"-V,(e, x) + a{-Xr (2.29) 

i=l 

where the independent functions Jj {i = 1 , . . . , n) are in involution and can be 
written as 

n 

Ji(e, x) = J2 Mi^)^] - (^(^ii^) with Ai{x) = g^{x)ai^{x). (2.30) 
In the case i = 1, we recover the Hamiltonian (j2.25|) . i.e., 

a=0 
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Proposition 2.4. The Jacobi system on T*£ defines a Stdckel system, with Stdckel 
matrix 

^^(^') = (-iry(^ (2-31) 

and potential functions 



^^^^^ = ^^-^^^ 



for i,k = 1, 



Proof. It is obvious from its definition that i? is a Stackel matrix. We just need to 
prove that A = B-^. To this aim we first prove a useful identity. Let us consider 
the integral in the complex plane 

1 [ u..w^^^ 



2i7r {z - A) U^{z) 

When R —>■ oo the previous integral vanishes because the integrand vanishes as 

for large R. We then compute this integral using the theorem of residues and we get 

the identity 

k=l ^ ' j^k 

Equipped with this identity let us now prove that 

n 
k=l 

Multiplying this relation by and summing over j from 1 to n, we get 

the equivalent relation 

n n 

k=i j=i 
which becomes, using ()2.30|) and ()2.27|) : 



k=l j^k 

Using the explicit form of g^{x) given in ()2.22j) and of the matrix B, this relation 
reduces to the identity ()2.33|) and this completes the derivation of ()2.3ip . 
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In order to get the functions as in ()2.1()|1 . let us resort to ()2.3()|1 and solve, 

for the unknown f], the following equation 



-aai{x) = '^Al{x)fj. 
Multiplying both sides by Bl, summing over i from 1 to n, and using ()2.3H) we get 

n n 

h = -a$:5^a.(x) = -^^5:(-ir(x'^r^-V.(a;) 

n 
i=0 



In view of flT^ . we have YJ'i=o{-^)\x^T~'(^i{x) = U^ii^'' - = 0, which 
completes the proof. □ 

Remark 2.5. 1. The fact that the geodesic flow on T*S is a Stackel system 
was flrst proved by Benenti in We have given here a new derivation, 
which makes the link between Moser's conserved quantities on and 
the Stackel conserved quantities on T*S. We have extended this link to the 
case where Jacobi's potential is admitted. 

2. Checking that the unconstrained observables Jj are in involution is most conve- 
niently done using their generating function (j2.26|) . Indeed it is easy to verify 
the relation 

{Gxi^, x), x)} = 0, A, /i G R, 

which implies, via ()2.29|) . and upon expansion in powers of A and /i, the rela- 
tions {Jj, Ij} = for any i, j = 1, . . . ,n. 

3. Some authors |21I22| have quantized the full set of commuting observables for 
the geodesic flow of the ellipsoid S C R"^^ in its unconstrained form, namely 
on T*R"+^. Notice though that in the reduction process from to T*S 
quantum corrections may prove necessary in order to insure self-adjointness of 
the quantized observables. Our point of view will be to perform the classical 
reduction in the flrst place and then to quantize the observables directly on 
T*S via a speciflc procedure that will be described in Section El 
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2.6 The Neumann system 

The Neumann Hamiltonian on {T*TV^~^^, Yl^=o d-Pa A dya) is 



^ = ^ E + «"^") (2.34) 



a=0 



with the real parameters < Cq < ai < . . . < a„. Under the symplectic reduction, 
with the second class constraints 

n n 

Zi{p,y) = Y^yl -1 = 0' ^2(p,2/) = Y.p-y'^ = O' (2.35) 

a=0 a=0 

it becomes a dynamical system on (T*S"', d^i A dx^). 

This system is classically integrable, with the following commuting first integrals 
of the Hamiltonian flow in T*R"^^: 



F.{p,y) = yl + Y. ^^''^' ^'^"^ ^ith a = 0,l,...,n. (2.36) 
The symplectic embedding 

i ' T* ' ^ j^n+l 

given by Zi{p,y) = and Z2{p,y) = preserves the previous conservation laws. 
Indeed the Poisson brackets of the restrictions Fa\rpt£ = Fa o l of the functions Fa 
are still given by the Dirac brackets (j2.18|) of the second class constraints (j2.35|) . 
This time we have 

n 

{Zi, z,} = -2j2yl^o, {Zi, Fa} = 0, 

which gives again 

{Fa\T*£, F^\rp,^} = 0. 

Let us introduce an adapted coordinate system on (T*S'", d^i/\dx^) much in the 
same manner as for the ellipsoid. 

We start with the following definition j^Hl of a coordinate system (x^, . . . ,a;") 

on 

yi _ nr=i(A-x^ 



- ^ na=o(^ - 
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The following inequalities hold: < ao < < ai < . . . < < a„. We get, in the 
same way as before, 



.2 _ ni=l('^a X* 



yi = \ (2.37) 

together with the following expression of the round metric g = Yla=o ^Vals" terms 
of the newly introduced coordinates, namely 

g = Y,^g,{x)idxr with g,(x) = -^|||| (2.38) 

with the notation ()2.2H) . Again, we put for convenience g*(a;) = l/gi(x). 

Our goal is to deduce from the knowledge of (j2.36|) the independent quantities in 
involution Ji, ...,/„ on {T*S^, d^i A dx*). The formula ()2.24|) relating unconstrained 
and constrained momenta still holds and yields the 

Proposition 2.6. The Neumann system {Fal-r* s"-) a=o n '^^^^^^'5 the following form 



F I 



where 

n n 
i=l j^i j=l 

Let us, again, posit 

n 

Gxit x) = Yl (-1)'-'A"-V,(e, x) + A" 
1=1 

where the independent functions Jj {i = l,...,n) are in involution and can be 
written as 

n 

m,x) = 5^ AKx)e| - cr,{x) with Ai{x) = g\x)al,{x), (2.39) 
i=i 

where the symmetric functions crj(x) are as in ()2.28|) . 
Using the relations 



(^i{x) = x\ and a^yl = - 

i=l a=0 a=0 i=l 

one can check that the Hamiltonian ()2.34j) is H = i/i. 
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Proposition 2.7. The Neumann flow on {T*S"',H) defines a Stdckel system, with 
Stdckel matrix 



k\n~ 



and potential functions 



for i,k = 1, . . . , n. 



Proof. To check that A = B~^, it is enough to use the identity (j2.33|) . The compu- 
tation of the potential functions fk proceeds along the same lines as in the proof of 
Proposition 12. 41 □ 

Remark 2.8. The involution property {/j, Ij} = for i,j = 1, . . . ,n, similarly to 
the case of the ellipsoid, is seen to follow from the relation {Ga(^, x), G^{^, x)} = 0. 

2.7 Test particles in generalized Kerr-Newman background 

Plebanski and Demianski have constructed in [211 EOl a class of metrics generalizing 
the Kerr-Newman solution in 4-dimensional spacetime. The former are also known 
as the Kerr-Newman- Taub-NUT-de Sitter solutions of the Einstein-Maxwell equa- 
tions. The metric, in the coordinate system {x^ , x"^ , , x^) = {p,q,a,T), retains the 
form 

^ idr + q^daf-^^idr-p^daf + ^4^dp^ + ^47^dq'' (2.41) 



p2 _|_ g2 p2 _|_ q,2 \ ^ ' ' X ^ ■ Y 

with 

X = 7 - ^2 ^ 2np - - - & y = 7 + - 2mq + eq^ - - g^ (2.42) 

3 3 

where (m, 7) are related to the mass and angular momentum of the Kerr black 
hole, (e, g) to the electric and magnetic charge; n is the NUT charge, and A the 
cosmological constant. The remaining parameter e can be scaled out to ±1 or 0. 

This metric, g, together with the electromagnetic field, locally given by F = dA 
where 



pz _l_ qz 



+ gp)dT + pq{gq — ep)da , (2.43) 
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provide an exact solution of the Einstein-Maxwell equations with cosmological con- 
stant A. Let us notice for further use that 

ViA' = 0. (2.44) 

Upon defining the 1-forms 



{dr — da) — \ ^ dq, 



M2 = J^^{dT + q^da), 



one constructs the 2-form 

y = pKAL-qMiA M2. (2.45) 

One can check that the twice-symmetric tensor P = —y"^, namely Pij = —yikyejs'^^ , 
is a Killing-Maxwell tensor (see ()2.7|) ). given by 

P = p^(K ®L + L®K) + q^{Mi ® Mi + M2 M2). (2.46) 

We thus recover Carter's result about the integrability of the Hamiltonian flow 
for a charged test particle in the generalized Kerr- Newman background in a different 
manner. 

Remark 2.9. The 2-form y in ()2.45|) defines what is usually called a Killing- Yano 
tensor pTllH]. 

The four conserved quantities in involution for the generalized Kerr-Newman 
system are, respectively, 

H = i g^^te - A,)(0 - A,), P = P^^(6 - A,)(0 - A,) (2.47) 

where P is as in ()2.46|) . and 

5 = 6-^3, f = U-A^. (2.48) 
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2.8 The Multi-Centre geodesic flow 

The class of Multi-Centre Euclidean metrics in 4 dimensions retain, in a local coor- 
dinate system (x*) = {t, {y"')) G R x R^, the form 

1 



V{y) 



{dt + AMdyy + V{y)^ (2.49) 



with 7 = 5ah dy^-dy^ the flat Euclidean metric in 3-space, and dV = ±-k{dA) where -k 
is the Hodge star for 7. These conditions insure that the metric ()2.49|) is Ricci-fiat. 

For some special potentials V{y), the geodesic flow is integrable as shown in 
[2IH UHl EH] • The four conserved quantities in involution are given by 

i/ = |g^^6e„ K = K%, L = L%, P = P'^^i^,, (2.50) 

where K and L are two commuting Killing vectors and P a Killing 2-tensor whose 
expressions can be found in the previous References. 

2.9 The Di Pirro system 

Di Pirro has proved (see, e.g., [2H]; P- 113) that the Hamiltonian on T*R^ 

2(7(x^, x"') + c(x'^)) 
admits one and only one additional first integral given by 

^ [c{x') {a{x\x')e, + h{x\ x')il) - 7(x\ x')^-^ . (2.52) 



(7(0;^, x^) + c(x'^)) 



In the case where the metric defined by H in (|2.51|) possesses a Killing vector, 
the system becomes integrable though not of Stackel type. This happens, e.g., if (i) 
c(x^) = const., or (ii) a = h and 7 depend on r = a/ {x^Y + (x^)^ only. 



3 A quantization scheme for integrable systems 

We wish to deal now with the quantum version of the preceding examples. Let us 
start with some preliminary considerations: 
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1. There is no universally accepted procedure of quantization, i.e., of a linear 
identification, Q, of a space of classical observables with some space of lin- 
ear symmetric operators on a Hilbert space. One — among many — of the 
pathways to construct such a quantization mapping has been to demand that 
the mapping Q be equivariant with respect to some Lie group of symplecto- 
morphisms of classical phase space. 

2. Similarly, there is no universally accepted notion of quantum integrability. 
However, given a classical integrable system Pi, . . . , P„ on a symplectic mani- 
fold {M.,iu), and a quantization mapping Q : Pj Pj, we will say that such a 
system is integrable in the quantum sense if [Pi, Pj] = for all i,j = l,...,n. 

3. A large number of integrable systems involve quadratic observables. We will 
thus choose to concentrate on this important — yet very special — case, both 
from the classical and quantum viewpoint. 

4. Among all possible quantization procedures, the search for integrability-pre- 
serving ones (if any) should be of fundamental importance. The quantization 
of quadratic observables we will present below might serve as a starting point 
for such a programme. 

3.1 Quantizing quadratic and cubic observables 

Let us recall that the space T\{M) of A-densities on M is defined as the space of 
sections of the complex line bundle |A"T*M|^ C. In the case where the configura- 
tion manifold is orientable, (M, vol), such a A-density can be, locally, cast into the 
form = /|vol|^ with / e C°°{M) which means that transforms under the action 
of a e Diff(M) according to / i— > a*/|(a*vol)/vol|'^. 

The completion H{M) of the space of compactly supported half-densities, A = i, 
is a Hilbert space canonically attached to M that will be used throughout this article. 
The scalar product of two half-densities reads 



where the bar stands for complex conjugation. 

We will assume that configuration space is endowed with a (pseudo-)Riemannian 
structure, (M, g); and denote by |volg| the corresponding density and by Ff^- the 
associated Christoffel symbols. 




19 



The quantization now introduced is a linear invertible mapping from the space 



of quadratic observables P = P2^{x)^j^k + + -^0(2;) to the space of second- 



order differential operators on H{M), viz A = P = (x)VjVfc + y4{(x)Vj + Ao(x)l 
where the covariant derivative of half-densities Vj0 = dj(f) — ^^^^(p (o^') locally, 
Vj0 = ((9j/)|volg|^) has been used. We furthermore require that the principal 
symbol be preserved (see below ()3.1|) . ()3.2|) and ()3.3|) ). and that P be formally 
self-adjoint, i.e., {(p/Pip) = (Pcpyip) for all compactly supported 0, G jFi(M). 
The quantization reads 

= -Pf (3.1) 
A{ = iP'i-WuPt (3.2) 



Ao = Po+'-VjP( (3.3) 



and admits the alternative form 



P = -V.oPf oVfc + ^ {P(oVj + VjoPi) + Pol (3.4) 
which makes clear the symmetry of the quantum operators. 

Remark 3.1. The formula ()3.4|) was originally used by Carter [TT] for proving the 
quantum integrability of the equations of motion of charged test particles in the 
Kerr-Newman solution. 

Remark 3.2. It is worth mentioning that formula ()3.4|1 actually corresponds at the 
same time to the projectively equivariant quantization [211 CS] and to the conformally 
equivariant quantization [13 Hn] Qo,i{P) '■ J^o{M) J^i{M) restricted to quadratic 
polynomials. 

One can check the relations: 

[Po, Qi] = t[Po, Qi]s = t{P^i}, (3.5) 



[Po, Q2] = -I (V,o[Po, Q2]i + [Po, Q2]'5°V,) = t{Po, Q2}, 



(3.6) 



[Pu Qi] = -\ (V,°[Pi, QiYs + [^1, Qi]5°V,) = ^{Pl, Qi}. 



(3.7) 
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Quantum corrections appear explicitly whenever k + £ > 2, as can be seen from 
the next commutators: 



where 



[Pl,Q2]=^{PuQ2} + ^Ap,,Q, 



(3i 



(3.9) 



is a scalar quantum correction that may vanish in some special instances, e.g., if the 
vector- field Pi is divergence-free (in particular if it is a Killing vector- field) . 

The previous formula can be found, in a different guise, in[TT]. Here, we will go 
one step further and compute the commutators [^2,^2] which involve third-order 
differential operators. To that end, we propose to quantize homogeneous cubic 
polynomials according to 

(3.10) 



P3 = ~2 °Va;°V£ + VjoVfcoP| oV< 

as a "minimal" choice to insure the symmetry of the resulting operator. 

Remark 3.3. The formula (|3.10|) precisely coincides with the projectively equi- 
variant quantization [7j Qo,i{P) '■ ^o{M) jFi(M) restricted to cubic polynomials. 

The previously mentioned commutator is actually given by 
[P2M 



[P2,Q2]?V,oVfcoV, 
+ \ (v,[P2,Q2]f) VfcoV, 



(3.11) 



+ 



V,Vfc[P2,Q' 



2J5 



+ 3 (VfeP. 



P2, 



where the skew-symmetric tensor 

p^y v.v^g'^i™ + p^[^p'i,^„,g"" - (p 



Q) 



(3.12) 



satisfies, in addition, Pp,Q = —Bq^p. We have used the following convention for the 
Riemann and Ricci tensors, viz R^'i^j^ = djVl^^ — {j ^ k) + . . ., and Rij = R^i^kj- 

We can rewrite the commutator ()3.11|) with the help of the quantization pre- 
scription ()3.4|) and ()3.10p as 



[P2,Q2]=l{P2,Q2} + lAp,,Q, 



(3.13) 
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where 

Ap,,Q, = -l{VkB'lQ,)^e (3.14) 

is a divergence-free vector-field associated with the tensor ()H.12|1 and providing 
the potential quantum correction for quadratic polynomials; recall that, according 
to (m, one has Ip^g, = {t/2){A'p^^Q^oVe + V,oA|,^ qJ. 
We thus have the 

Proposition 3.4. The commutator of the quantum operators P and Q associated 
with two general quadratic polynomials P = P2 + Pi + Pq and Q = Q2 + Qi + Qo 
reads 

-[P, Q] = {P^} + Ap,,Q, + lp,,Q, - Aq,,p, (3.15) 
where the third-order differential operator {P, Q} is given by i3.1(J\} . 

Proof. The formula ()3.15p results trivially from the previously computed commuta- 
tors and from collecting the anomalous terms appearing in (j3.8|) and (j3.13|) only. □ 

Remark 3.5. In the special case where Q2 = H a.s given by ()2.3|1 . the anomalous 
tensor ()3.12|) takes the form 

and reduces to 

= -P'imf (3.16) 

if P is a Killing tensor ^T] . 

Remark 3.6. In the particular case where H = ^g^^{C,j — eAj){C,k — eAk) is the 
Hamiltonian of the electromagnetic coupling, our quantum commutator (|3.15|) re- 
duces to Carter's formula (6.16) in [TT| . 

The purpose of our article is, indeed, to study, using explicit examples, how clas- 
sical integrability behaves under the "minimal" quantization rules proposed in [TT] 
and somewhat extended here. The next section will be devoted to the computation 
of the quantum corrections in ()3.8|) and ()3.13|) for all the examples that have been 
previously introduced. 
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3.2 The equivariance Lie algebra 



So far, the transformation property of the quantization rules ()3.4p and ()3.10|) under 
a change of coordinates has been put aside. It is mandatory to investigate if these 
rules are consistent with the map Q : P t-^ P (which has been defined for cubic 
polynomials, P = X]fe=o ''^n • • • '^«fc' '^^^j) being equivariant with respect to 
some Lie subgroup of the group of diffeomorphisms of configuration space, M. 

Restricting considerations to the infinitesimal version of the sought equivariance, 
we will therefore look for the set q of all vector fields X with respect to which our 
quantization is equivariant, namely LxQ = 0. From its very definition, g is a Lie 
subalgebra of the Lie algebra, Vect(M), of vector fields of M. The previous condition 
means that, for each polynomial P, the following holds: 

Lx{Q{P)<l>) - Q{LxP)<P - Q{P)Lx<P = (3.17) 

where Lx0 denotes the Lie derivative of the half-density cf) oi M with respect to the 
vector field X G and LxP = {X, P} is the Poisson bracket of X = and P. 

Let us recall that, putting locally = /|vol|5 g JFi with / G C°°(M), we get 
the following expression for the Lie derivative: Lxcj) = {Xf + idiv(X)/)|vol|^, or 
with a slight abuse of notation, Lx(p = X^Vj4> + i(VjXJ)0 = ^{X^oVj + VjoX^)0, 
that is 

Lx(j) = -X(j) (3.18) 
i 

for any X G Vect(M). 

The equivariance condition ()3.17p must hold for any G JFi and thus translates 

into 

[X,P]=t{XJ^} (3.19) 

for any X E Q and any cubic polynomial P. The Condition (j3.19j) characterizes 
the Lie algebra g we are looking for. We will consider successively the case of 
polynomials of increasing degree: 

(i) Returning to the previous relations (j3.5|) . (j3.7p together with X = Pi and 
P = Qo + Qi, we readily find that the Lie algebra gi spanned by the solutions 
of (j3.19j) restricted to polynomials P of degree one is gi = Vect(M). 

(ii) Let us now proceed to the case of quadratic polynomials P = P^^^j^k- The 
relations (|3.7p and (j3.9p give, in that case, the following equivariance defect 

[X, P] - 100"} = '-VjoP^KVk{ViX')l. (3.20) 
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This defect vanishes for any such P iff VkCVeX^) = 0, i.e., 

d{dw{X)) = 0. (3.21) 

The vector fields X with constant divergence span now a subspace 92 C gi which is, 
indeed, an infinite dimensional Lie subalgebra of Vect(M). The "minimal" quanti- 
zation restricted to quadratic polynomials is therefore equivariant with respect to 
the group of all diffeomorphisms which preserve the volume up to a multiplicative 
nonzero constant. 

(iii) Let us finally consider homogeneous cubic polynomials P = P-'^^^jCk^e and 
compute the equivariance defect in this case. A tedious calculation leads to 

[X, P] - i0O^} = iZ, Z = Z^ij, (3.22) 

with 



Z^ 

where 



P^'^V^div(X) - P^™[^LxrJJ„ (3.23) 



LxLL = V.V^X' - R',X- (3.24) 



is the Lie derivative of the symmetric linear connection V with respect to the vector 
field X. 

Proposition 3.7. The Lie algebra g C Vect(M) with respect to which the "minimal" 
quantization and \3.1U\} is equivariant is aff(M, V), the Lie algebra of affine 

vector fields of (M, V) . 

Proof. The equivariance condition ()3.19p , defining the Lie algebra 03 we are looking 
for, is equivalent to Z = in (j3.22p for all symmetric tensor fields P-?^^, i.e., thanks 
to (IT^ to 

Tf V,div(X) - T'-H^vt = 
for all tensor fields T^"^-' = Tlf'^''\ This readily implies that 

26l6^Vm)div{X) + 5;X;,rL) - ^Ji^xFL) = 0. 

Summing over i = j, one gets 

2n5^V,div(X) + 45,^V^div(X) + (n + 1)LxTL ' €LxT\, - S'^LxT^, = 0, 

where n = dim(M), hence Vjdiv(X) = and Lx^fj = S^(fj + 6j(fi for some 1-form (p 
depending upon the (projective) vector field X. The expression ()3.24p of the Lie 
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derivative of the symmetric connection V then yields Lx^jj = {n + l)^Pi = since 
we have found that ViVjX^ = 0. This entails LxT^j = 0, proving that g = 03 is 
nothing but the Lie algebra aff (M, V) of affine vector fields. □ 

We thus obtain the nested equivariance Lie algebras 

= aff (M, V) C 02 C 01 = Vect(M) 

where 02 is the Lie algebra of vector fields with constant divergence. (Note that 
if M is compact without boundary, 02 reduces to the Lie algebra of divergence-free 
vector fields.) 

Conspicuously, our quantization scheme turns out to be equivariant with respect 
to a rather small Lie subgroup of Diff (M), namely of the affine group of (M, V). It 
would be interesting to investigate to what extent the equivariance under the sole 
affine group, GL(n, R) ix R", of a flat affine structure (M, V) allows one to uniquely 
extend to the whole algebra of polynomials the quantization scheme we have devised 
for cubic polynomials. 

3.3 The quantum Stackel system 

The quantization of the general Stackel system (see Section lT^ has first been under- 
taken by Benenti, Chanu and Rastelli in We will derive, here, the covariant 
expression of the quantum correction associated to the "minimal" quantization, with 
the help of the results obtained in Section 13.11 

Denote by = /2,i -|- Io,i the z-th Stackel conserved quantity, i = l,...,n, 
in 1)2.111) where the indices and 2 refer to the degree of homogeneity with respect 
to the coordinates ^. Applying (|3.15p with Pi = Qi = 0, P2 = l2,i and Q2 = hj 
one gets 

Remark 3.8. This result shows that there are no quantum corrections produced 
by the potential term. More generally, start with a system defined by independent, 
homogeneous, quadratic observables Hi, ... , Hn which is integrable at the classical 
and quantum levels. Consider a new set of observables Hi + Ui, . . . , Hn + Un obtained 
by adding potential terms Ui, . . . ,Un; if the new system is classically integrable, it 
will remain integrable at the quantum level. 

We are now in position to prove the following 
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Proposition 3.9. The quantum correction hS.l^) of a general Stdckel system, with 
commuting conserved quantities Ji, . . . , J„ defined by \2.11\) . retains the form 

B^l,^^ = -2lf^RstlZ (3.25) 

for i,j = l,...,n, where Rgt denotes the components of the Ricci tensor of the 
metric associated with the Hamiltonian Ii . 

Proof. As a preliminary remark, let us observe that the Stackel metric, given by 
()2.1U|) . needs not be Riemannian. So we will write it 



i\2 



t'-^r^'Ue-T- (3.26) 

i=l ' 



a=l 



where {9°' = dx"" / \/\A^\)a=i,...,n is the orthonormal moving coframe and the signature 
of g is given by rja = sign(A"). We will denote by {ca = ^/\Af\da)a=i,...,n the 
associated orthonormal frame with respect to the metric rjab = rjaSab used to raise 
and lower frame indices. 

Let us recall, in order to fix the notation, that the connection form u satisfies 
the structure equation d9°- + cj"^ A 9^ = and the associated curvature form, Q, 
given by = duj°'fj + u"'^ A tu'^^, is expressed in terms of the Riemann tensor by 
= IR\ 9"" A 9'^. The indices a,...,d run from 1 to n and the Einstein sum- 
mation convention is used when no ambiguity arises. Denoting by R^i^k the local 
components of the Riemann tensor, we have R°'i,^cd — ^^ijk^l^c^d- 
We start off with the calculation of the connection form, u, and of some compo- 
nents of the curvature form, Q. Straightforward computation, using relation (|2.13|) . 
then yields for the non-vanishing components of the connection 

1^1 1 

the other nontrivial components ujab,b are obtained accordingly. For the curvature, 
a lengthy computation gives the special components 

Rac,cb = 3 { — riaWca,ci^ab,a — Vb^cb,c^ba,b + Vc^^ca,c^cb,c) , a ^ b, (3.27) 

which will be needed in the sequel. 

Two last ingredients are the introduction of the frame components of various 
objects. We will denote the Killing tensor J2,i (resp. I2j) as P (resp. Q). Their 
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frame components P = P Cb <S) Cc, and similarly for Q, will be 



P^"" = PbSbc, Pb 



Q''" = QbSbc, Qb 



A'] 



(3.28) 



2\A\y ^ ^""'^ 2\A\[ 
The covariant derivative will have the frame components 

T^cPab = ^c{Pab) — ^^a,cPsb — ^^hfiPas- 

The equations which express that P"^ is a Killing tensor are now 

eb{Pa) = '2Uab,aiVaPa - VbPb) , a ^ b, 
ea{Pa) = 0, 

where the repeated indices are not summed over. One can check that they hold true 
using the explicit form of Pa given in (|3.28|) and the identity (|2.13|) . 

Using all of the previous information one can compute the frame components 
of the various pieces appearing in the tensor Bpg. We have successively 

P^['V.VtQJ]* - (P ^ Q) = 



(3.29) 



PiQj - ViPiViQj + ViQiViPj - {i^ j] 



and 



PiQj - ViPiViQj + miViPj - (« ^ j] 



Combining these relations, and using ()3.27|) . we get 

P^f'V.VtQ^l* - (P ^ g) - V.P*['VtQ^> = 



PiQj - viPiViij + vmviPj - ^ J) 



ViPiViQj - ViQiViPj - 



Let us then compute 

P^t'P^UQ™ - Q) = i^r^,P.M. 

I 

Collecting all the pieces leaves us with 

P^l^V.VtQ-^'l* + P^l^P^'U.Q™ - (P ^ g) - V.P*[*VtQ^> 

hJ2iViRmjiPiQj -PjQi)- 
27 



(3.30) 



The last sum is nothing but the frame components of the tensor —P'^^^RstQ-'^^, so 
that we have obtained the tensorial relation 

P'^'VsVtQ'^' + P^['i?^l^,,Q™ - (P ^ g) - VsP'^'VtQ'^' = -P'^'RstQ'^', (3.31) 

which implies 

= -2P'^'RstQ'^\ (3.32) 

in agreement with j3] . This ends the proof of Proposition 13.91 □ 

Now we can come to the central point of our analysis: is a Stackel system 
integrable at the quantum level? The answer is given by the following 

Corollary 3.10. ('Q l^j A Stackel system is integrable at the quantum level iff 

Rij = for i 7^ j, where i, j = 1, . . . ,n, (3.33) 

in the special coordinates which are constituent to this system. 

Proof. The Killing tensors J2,j are diagonal, for i = 1, . . . , n, in the Stackel coordinate 
system, and the proof follows from ()3.25|) . □ 

The conditions ()3.33|1 are known as the Robertson conditions [HI] , as interpreted 
by Eisenhart ^H]- Quite recently, Benenti et al 0] have refined the definition of 
the separability of the Schrodinger equation and shown that, for Stackel systems, 
the Robertson conditions are necessary and sufficient for the separability of the 
Schrodinger equation. As mentioned in Remark 12.11 the classical integrability is 
equivalent to the separability of the Hamilton- Jacobi equation; the situation for 
these systems can be therefore summarized by the following diagram: 

Classical integrability <^==4> separable Hamilton- Jacobi 
JJ. provided Rij = {i ^ j) 
Quantum integrability <^=^ separable Schrodinger 
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3.4 The quantum ellipsoid and Neumann systems 

It is now easy to prove that the elhpsoid geodesic flow (see section including 
the potential given in ()2.1fi|l . is integrable at the quantum level. Using the coor- 
dinates (x*) and the (Riemannian) metric given by ()2.22|) . one can check that the 
Ricci tensor has components 

_ A/^ >^ J_ apai ■ ■ - an 

~ ^ x^- ■■X"' ' 

and therefore satisfies the Robertson conditions. As already emphasized, the oc- 
currence of an additional potential is irrelevant for the quantum analysis since the 
potential terms do not generate quantum corrections (see Remark l3.8|) . 

Similarly we get the quantum integrability for the Neumann system (see Sec- 
tion EIEI) using the metric on given by ()2.H8|1 . The Ricci tensor being given 

by 

Rtj = {n- l)gij, 
the Robertson conditions are again satisfied. 

3.5 The quantum generalized Kerr-Newman system 

The quantization of the four commuting observables ()2.47|) and (j2.48p is straight- 
forward. 

In view of the relations given in Section El all quantum commutators vanish 
except for [i?, -P]; this is due to the fact that the conserved quantities S and T 
(see (I2.48p ) are Killing- Maxwell vector fields. 

The anomalous terms in the previous commutator are Ap^^H2^ ^Pi,H2 ^"^^ ^P2,Hi 
where P2 = P^^'^^,, ^2 = lg''^i^„ Pi = -2P''^^A, and H, = -g'^^^A^. 

The vector field Ap^ u^ given by ()3.14p actually vanishes because, cf. ()3.16p . 
B^p^ _H"2 ~ —P^^-'P^i = as a consequence of 1)2.7^ : indeed the tensor P anti-commutes 
with the electromagnetic field strength F, implying that it commutes with the stress- 
energy electromagnetic tensor, hence with the Ricci tensor in view of the Einstein- 
Maxwell equations [TT] . 

The two other anomalous terms ()3.9p also vanish as it turns out that VjA-^ = 
(see (ESI)) and Vj{P^^Ak) = 0. 

This derivation reproduces and extends Carter's results to the generalized Kerr- 
Newman solution, in a somewhat shorter manner. 
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Remark 3.11. Our analysis of quantum integrability for the generalized Kerr- 
Newman solution in 4 dimensions can be carried over into recent work 1221 
dealing with 5-dimensional black holes. In these cases, classical integrability follows 
from the existence of 3 Killing vectors and 1 quadratic Killing tensor, besides the 
Hamiltonian. These metrics being Einstein, the above arguments given for the 
generalized Kerr-Newman case apply just as well, insuring quantum integrability. 
This fact is in agreement with the separability of the Laplace operator. 

3.6 The quantum Multi-Centre system 

For this example too, the quantization is straightforward. The single point to be 
checked for quantum integrability is just the commutator [H,P], with the possible 
quantum correction ()3.1(ij) given by —P^^^uf. Here it vanishes trivially since these 
metrics are Ricci-flat. 

3.7 The quantum Di Pirro system 

As seen in Section 12. 9[ the classical integrability of this system is provided by three 
commuting observables: on the one hand H, P respectively given by ()2.51|) and 
(I2.52p . and T = ^3 if c(x^) = const., and on the other hand H, P and J = C,ix'^ — C,2X^ 
if a = 6, 7 depend on r only. 

At the quantum level, the Killing vectors T and J do commute with H according 
to ()3.8|) and ()3.9|) . As for the commutator [P,H] of the quantized Killing tensors, 
it is given by (jSISl), namely Bp^n = -^P^^^Pf dj A dk, and one finds 

16 (7(x^, x^j + c(x'^))'^ 
+ h{x^,x^)d2i{x^,x^)d2Kd-^. 

For the system {H,P,T), this quantum correction vanishes since c'(x^) = 0, 
implying quantum integrability. However, for the system {H,P,J), in the generic 
case 7 7^ const., we get Bpjj 7^ 0, showing that the minimal quantization rules may 
produce quantum corrections. 
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4 Discussion and outlook 



It would be worthwhile to get insight into the status of our "minimal" quantization 
rules and to their relationship with other bona fide quantization procedures. Among 
the latter, let us mention those obtained by geometric means, and more specifically 
by imposing equivariance of the quantization mapping, Q, with respect to some 
symmetry group, G, e.g., a group of automorphisms of a certain geometric structure 
on configuration space, M. We refer to the articles [211 HSJ Uni HZl E] for a detailed 
account on equivariant quantization. The two main examples are respectively the 
projectively, G = SL{n + 1, R), and conformally, G = 0{p + 1, g + 1), equivariant 
quantizations which have been shown to be uniquely determined fI^ [T7\ [W\ [W\ . 
For instance, the conformally equivariant quantization Qi : !Fi{M) J-'i{M) has 
been explicitly computed for quadratic [12] and cubic [2S1 observables; for example, 
if P = P^^^i^j we then have 

Qi(P) =P + (3:,V,V,{P'^)+l3^g'^gutV,V,{P''')+l3,R,,P'^ +l3^Rg,,P'^ (4.1) 

where the "minimal" quantum operator 

P = -S/ioP'KVj (4.2) 

is given by (Hnil . together with /?3 = -n/(4(n + 1)), (3^ = -n/(4(n + l)(n + 2)), 
/?5 = nV(4(n-2)(n + l)), = -n^ /{2{n^ - A){n^ -1)), assuming n = dim(M) > 2. 
In M.lj) we denote by Rij the components of the Ricci tensor and by R the scalar 
curvature. The formula ()4.H) provides a justification of the term "minimal" for the 
mapping P P given by (|3.4|) and (j3.1U|) . 

We have checked that, in the special instance of the geodesic fiow of the ellipsoid 
discussed in Section 12.51 the quantum commutators of the observables /j defined 
in ()2.3m) . namely [Qi (li), Qi (/j)], fail to vanish for z 7^ j = 1, . . . , n. Had we started 
from the expression 1)4.11) with adjustable coefficients P3, . . . , Pq, the requirement that 
the latter commutator be vanishing imposes /Ss = . . . = /Sg = 0, leading us back to 
the minimal quantization rule ()4.2|) . 

Despite their nice property of preserving, to a large extent, integrability (from 
classical to quantum), the "minimal" quantization rules still remain an ad hoc pro- 
cedure, defined for observables at most cubic in momenta, and do not follow from 
any sound constructive principle, be it of a geometric or an algebraic nature. The 
quest for a construct leading unambiguously to a genuine "minimal" quantization 
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procedure remains an interesting challenge. As discussed in Section 13 .21 the equiv- 
ariance assumption with respect to the affine group might be helpful for determining 
the sought "minimal" quantization of polynomials of higher degree. This analysis is 
required for the quantization of, e.g., the newly discovered integrable systems [Hj 
which involve cubic Killing tensors. 

Another field of applications of the present work could be the search for quantum 
integrability of the geodesic flow on the higher dimensional generalizations of the 
Kerr metric which have been lately under intense study fl9 \ IT2 | l36j. 

Still another perspective for future work would be to generalise the previous 
computation of quantum corrections to the case of classical integrability in the 
presence of an electromagnetic field in a purely gauge invariant manner. In particular 
the approach presented in Section 12.21 should be further extended at the quantum 
level via the quantization of the Schouten-Maxwell brackets. 
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